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Abstract. We study the partition-function zeros in mean-field spin-glass models. 
We show that the replica method is useful to find the locations of zeros in a complex 
parameter plane. For the random energy model, we obtain the phase diagram in 
the plane and find that there are two types of distribution of zeros: two-dimensional 
distribution within a phase and one-dimensional one on a phase boundary. Phases with 
a two-dimensional distribution are characterized by a novel order parameter defined 
in the present replica analysis. We also discuss possible patterns of distributions by 
studying several systems. 
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1. Introduction 

Spin-glass (SG) systems are known to sliow exotic thermodynamic properties and have 
been studied for many years [HOIS]. The existence of many quasi-stationary states leads 
to a nontrivial critical behavior at low temperatures. Other systems such as structural 
glasses and other theories like information processing share these properties. Therefore, 
a more general understanding is required for this randomness-induced transition. 

As a general method to understand the critical phenomena of statistical mechanical 
systems, Lee and Yang proposed to study the zeros of the partition function jH |5] . Since 
the partition function is positive by definition, the zeros are located in the unphysical 
region of the complex parameter plane. In finite systems, zeros are far from the real 
axis. However, as we make the system size large, they reach a point on the real axis 
at the thermodynamic limit. If the partition function goes to zero at a certain point 
in parameter space, the corresponding thermodynamic functions are singular, and this 
point is identified as the phase transition point. For example, in the pure ferromagnetic 
Ising model, the circle theorem states that the zeros form a unit circle on the complex 
fugacity plane [5^, ^ . In the same way, we can consider zeros of the partition function 
not only for complex fields but also for other complex parameters. For example, the 
zeros in the complex temperature plane are known as the Fisher zeros [7]. 

While there are many studies on the partition-function zeros for pure systems, 
relatively few results exist for random ones. The zeros of the ±J Ising model were 
examined numerically in [HI |9l [101 E] for finite-dimensional systems and in [12j on the 
Bethe lattice. The results show that the zeros are distributed in the complex plane 
in a very different manner compared to pure systems. That is, they are distributed 
two dimensionally rather than one dimensionally. The two-dimensional distribution 
approaches the real the system size is made large, which is considered to be the 

onset of the SG phase transition. Therefore, it is important to understand how such a 
distribution is formed in the SG systems. However, most of the previous studies were 
done numerically, and we need a reliable analytical method. 

As a first step to understand the distribution of zeros in SG systems analytically, 
we propose a new method to use replicas for mean-field SG models. The theory of 
spin glasses has been well understood in terms of mean-field models with infinite-range 
interaction. It was found that the replica method is useful for analytic evaluation 
of the average free energy, and replica symmetry breaking (RSB) corresponds to the 
SG state [131 [IH IIS]- The advantage of studying mean-field models is that analytical 
calculations are possible, and we can closely study how such a SG state is obtained. 
Therefore, it would be useful to study the zeros of the partition function in mean-field 
SG models. 

As a matter of fact, the analytical result is already available for the simplest SG 
model, the random energy model (REM) [161 113 IIH]- In this model, the distribution 
of zeros in the complex temperature plane was calculated analytically in [19] and 
confirmed numerically in [20] . The case of complex magnetic field was also studied both 
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analytically |21j and numerically [22]. The method used in [T^ [21] was to calculate 
the density of states for a given energy and is very specific to the REM. In fact, the 
spin degrees of freedom have not been treated, and it is hard to extend the calculation 
to other spin systems. Therefore, we propose a more general and systematic method 
using replicas so that the extension to other models is possible in principle. The use 
of the replica method has great advantages since many technical methods have been 
estabhshed, and we can utilize various properties obtained in the previous studies. 

The organization of this paper is as follows. In section [21 we give a brief review 
of partition-function zeros. Throughout this paper, we treat the REM and its variants. 
The model and the method of calculation are described in section [3l Then, in section [H 
we consider the Fisher zeros of the REM and re-derive the result obtained in [19]. Most 
of our main ideas and the essence of the calculations are contained in this section. The 
following sections are the applications of the method to several models. We treat the 
generalized REM (GREM) [231 [22 [23 [21] in section [5]to study the effect of the RSB for 
the distribution of zeros. Also, to study the first-order transition, we consider the REM 
with ferromagnetic interaction in section [6] The system in magnetic fields is considered 
to find the Lee- Yang zeros in [71 Finally, in the last section [HI we give conclusions and 
discuss issues for further study. 



2. Partition-function zeros 



We give a brief introduction of the partition-function zeros to derive several formulae 
used in the following sections. We treat Ising spin systems in a magnetic field h and the 
Hamiltonian is generally written as 

N 

H = Ho-hY,S„ (1) 

i=l 

where Si is the spin variable on site i, N the number of spins and Hq is the field- 
independent part of the Hamiltonian. The partition function can be expressed by using 

y = e^^^ as 

Z = Tre-'^^»+^'^i:"i^» = e-""^^ [a, + a,y + a^y^ + ■ ■ ■ + a^y^) , (2) 

where (5 is the inverse temperature and oq, ai, ■ ■ ■ , are coefficients determined by Hq. 
Apart from the overall factor e~^^'^, this function is a polynomial of the iVth degree in 
y and is factorized as 

N 

Z = a,,e-''^'^Y[{y-y{z)), (3) 

i=l 

where y{i) are the complex numbers and represent the zeros of the partition function. 
Then, the free energy per spin / = — (l/A^/3) InZ can be written as 

- /3f = -^InaN - Ph + J dzidz2 py{zi, Z2)\n{y - z), (4) 
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where yi and y2 denote the real and imaginary parts of y, respectively, and the same is 
for zi and Z2- The density of zeros Py{yi,y2) is introduced as 

1 ^ 

1=1 

Thus, the free energy can be written as an integral over zeros in the complex y plane. 

The expression of the free energy in terms of the density of zeros indicates the 
following features on the analyticity of thermodynamic functions. If the free energy (jlj) 
has a singularity, it should come from the integral at the point z = y. Since y = e^^^ is 
real and zeros are not on the real axis for finite systems, the phase transition indicated 
by the singularity is realized only at the thermodynamic limit. We can observe how 
zeros approach the real axis by increasing the system size A^. 

Our analysis in the following sections is based on the formula 

+ (6) 

which is derived from the expression (|3]) [19]. Here, the magnetic field in the partition 
function takes a complex value. That is, the density of zeros can be obtained from the 
absolute value of the partition function with complex-valued fields. 

We mainly consider in the following the Fisher zeros, the partition-function zeros 
in the complex temperature plane. Assuming that the partition function is factorized 
with respect to (3 as 

Z{P)=CX{{(5-m. (7) 



we can write 



hiZ{(5) = \YiC + N j dzidz2Pisizi,Z2)\n{l3- z), (8) 
/32) = TTl^ + In \Z{(3 = (3, + ^Z^^)!, (9) 



27riV \dl3l dp, 



where is the density of zeros in the complex-/? plane. We note that pp is not 
normalized to unity. This is because the partition function is a polynomial of infinite 
degree in /3 even for finite A^. 

3. Mean-field SG models and replica method 

3.1. Model 

We treat the REM which is represented by randomly distributed energy levels. For a 

given set of energies E2 - ■ ■ , E2N), the partition function is written as 

2N 

Z = Y^e-^'^\ (10) 

i=l 

Each energy level is taken from the Gaussian distribution 
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It is well known that this model is equivalent to the p — t- oo limit of the p-body 
interacting Ising spin model whose Hamiltonian is given by 

N 

H = — ^ Ji-^i^...i^Si-^Si^ ■ ■ ■ Sip, (12) 

il<i2<-<ip 

where the probability distribution of the interaction Ji^-.-i^ is Gaussian as 



Pj('^«i-v) = \/^exp|--j^(j,,,,,...g |. (13) 

The REM is solved exactly, and we can find a phase transition between the 
paramagnetic (P) and SG phases at /3 = /3c where 



/3cJ = 2Vln2. (14) 

At temperatures lower than Tc = l//3c, the system freezes to its ground state and the 
entropy goes to zero. This SG state is obtained by the one-step RSB (IRSB) ansatz in 
the replica formalism. 

3.2. Replica method 

In order to treat disordered systems, it is useful to take an ensemble average over disorder 
parameters. The self-averaging property ensures that one can study the thermodynamic 
limit of the system by using the averaged free energy. Then, the replica method is useful 
to calculate the average. As we described in the previous section, in order to obtain 
the density of zeros, we need to calculate [In |Z|] where the square bracket denotes the 
average. The main idea of this paper is to use the formula 

ln|Z| = lim^ . (15) 

Compared to the standard procedure to treat [Z^] for the free energy —f3F = In Z, we 
see that the number of replicas is doubled: = Then, the resultant SG 

order parameter Qab is represented by a 2n x 2n matrix. With complex parameters, the 
elements of the order-parameter matrix can take a complex value. However, in the case 
of the REM, they are real, which allows us to obtain the result easily. 

Since the replica space is doubled and the parameter is complex, we can find new 
phases which do not exist in systems with real parameter. In such a phase, it is shown 
in the following section that the factorization property 

P|2n] = (16) 

does not hold. In fact, this property is closely related to whether the two-dimensional 
distribution of zeros appears or not. 

4. Fisher zeros of the REM 



As a first example, we consider the REM with complex-/?. In this case, the density of 
zeros has been obtained analytically by Derrida in 1991 [19]. His method corresponds 



Replica analysis of partition-function zeros in SG models 



6 



to the microcanonical calculation in some sense since the density of states as a function 
of energy is calculated to obtain the canonical partition function. Here, we re-derive the 
result by the canonical replica calculation. 

4.I. Replica method 

For a given [Ei, E2, - ■ ■ , E2N), the partition function is defined as ( ITOl) . Then, we can 
write 

1^1'^ = E E^^p -PT.^AM)E, - p*Y.^,{K})E, , (17) 

KIR} I i=i i=i J 

where 

n 

'^.({U) = E^W' (18) 

a=l 

and the configuration sums are taken over {ia} = {ii, «2, ■ ' ' ; "^n) and {i'^} = {i[, i'^-, ■ ■ ■ ,in) 
with each ia and i'^ taking 1, 2, ■ ■ ■ , 2^. Taking the average ( 1TT|) . we obtain 

(19) 

Here, three order parameters are introduced as 

Qab = (20) 

q'a, = h^i^ (21) 

Qab = Si^i'^, (22) 
which can be written in a 2n x 2n matrix form as 

Q-{y^.)- (23) 

The n X n matrix q represents the overlap between configurations in {ia}, q' in {i'a} 
and q in {ia} and {i^}- Each matrix element takes or 1. At the thermodynamic limit 
N oo, the saddle-point configurations dominate the sum, and we can write 

I \ ah ah ab / J 

(24) 

where Nsq is the entropy function for configurations to give the saddle-point solution. 
Since is real, the saddle-point solution must satisfy the relation 

E?-^' = E^1^- (25) 

ab ab 



Then, 



\Zn = exp jiV i^s, + y (/3? - /3l) X; qab + y (/?? + /3,^) E ^ab^ } • (26) 
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As possible saddle-point solutions, we can consider the replica symmetric (RS) and 
IRSB ones. Among them, we choose the solution that gives the maximum ln|Z|. We 
also impose the condition that the thermodynamic entropy must be nonnegative. The 
definition of the entropy for complex-/? is nontrivial, and we discuss it in the next 
subsection. 



4-.2. Entropy for complex temperature 

In the standard theory of statistical mechanics, the partition function is written as 

Z = JdEe-^'^+'^^l (27) 

The entropy S{E) is defined as the logarithm of the number of states for a given energy 
E. At the thermodynamic limit, the saddle point (3 = dS/dE contributes to the energy 
integral. When {3 is complex as /3 = /3i + i/32, the saddle-point energy is also given by 
a complex number. We change the integral path so that the path crosses the saddle 
point. For an obtained saddle point E = Ei + iE2, we can write 

|^|2 ^ ^-2l3iEi+2l32E2+2SiEi,E2) (^28) 

where we introduce the real function S{Ei, E2). We write for $ = In |Z|, 

132) = -PiEi + /32^2 + S{E,, E2). (29) 

Thus, <l'(/3i,/32) is obtained from S{Ei,E2) by the Legendre transformation, and vice 
versa. The parameters /3i_2 and Ei^2 are related as 

^' - -IK' ^= ^ W. <^»' 

S{Ei, E2) defined here is considered a natural generalization of the ordinary entropy 
for real parameters. It is real and we employ the nonnegative condition of S{Ei,E2) 
when we determine the effective domain of a phase. 



4-3. Saddle-point solutions and phase diagram 

Now we consider the saddle-point solutions for fl26l) . Since qaa = I'aa — ^ ^y definition, 
we see that the simplest solution is given by 

Qab = q'ab = Sab- (31) 

This result is known as the ordinary RS solution and corresponds to the P phase. 
However, (I3T1) does not completely specify the saddle point, and we have an additional 
order-parameter matrix q which represents overlaps between the configurations {ia} and 
{i'a}- If there is no correlation between them, we have 

qab = 0, (32) 

which we call the PI phase. It is also possible to consider 

qab = Sab, (33) 
I More generally, this should be written as ^^j^ (jab = n since depends only on the sum. 
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which means that the configurations are completely the same: {ia} = {i'a}- We call 
the corresponding phase the P2 one. These solutions can be graphically expressed in 
figure [H n states are chosen from 2^-configurations. Figure [U^a) represents the PI 
phase and [U^b) the P2 one. In principle, we can consider intermediate states between 
PI and P2 phases such that < J2ab Qab < n, but they do not give the maximum In |Z|. 

In the PI and P2 phases, the numbers of states giving the saddle-point value of Q 
are respectively given by 

{„2Afnln2 pi 

We have for = [ln|Z|]/iV 



ln2 + i(/32-/3|)j2 PI 



^ |ln2 + i/32j2 P2 • 

The entropy in each phase is calculated from (1291) . We obtain the entropy per spin 

r ln2-i(/3?-/3|)J^ PI .3g. 

The condition s > gives us 

-Pf + Pl + P^^>0 (37) 
in the PI phase and 

|/3i|</3c/2 (38) 
in the P2 phase. Here, the inverse critical temperature /3c is defined as ([1 



(35) 



Next, we treat the IRSB SG solution which can be expressed in figure [2l We have 
for Qab and q'^f, 

II 9a6 = ^ q'ab = mn, (39) 

ab ab 

where m is the Parisi breaking parameter. For the configuration in figure [2]^a), there is 
no overlap between {ia} and {i'^} and 

qab = 0. (40) 
For that in figure [2](b), (jab has the IRSB structure and 

ab 

Then, for each saddle point. 




In2 + ^(^2_^2)j2 



ln2 I m o2 72 • ) 



(43) 
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Figure 1. RS saddle-point config- 
urations pip . TliG liorizontal axis 
represents the index of configura- 
tions running from 1 to 2^. The 
vertical axis represents the num- 
ber of each configurations, (a) PI 
phase jSH) and (b) P2 ([33|. 
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Figure 2. IRSB saddle-point 
configurations (|39|) . (a) SG phase 
with q = and (b) that with q ^ 0. 



and we obtain 

(44) 




2 



In both the cases, the entropy is shown to be zero. The former solution is always smaller 
than the latter one and is discarded. 
We summarize the result as 

^{Pl-Pl + Pl) {Pi - PI < PI) PI 

f(/3? + f) (|/3i|</3c/2) P2 . (45) 

^Pc\Pi\ SG 

We compare these solutions to take the maximum value of 0. The condition that PI is 
larger than P2 is 

Pl + Pl<^. (46) 
and that PI is larger than SG is 

m-Pcf>pi (47) 

SG is always larger than P2. They are equal when \Pi \ = Pc/2. In conclusion, the phase 
diagram is obtained as in figure [31 
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Figure 3. Phase diagram and 
distribution of partition-function 
zeros of the REM in the complex-/? 
plane. Zeros are distributed in the 
shaded area and on the bold lines. 
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Figure 4. The specific heat 
calculated from ([53l) . The solid 
line represents the exact result c = 
{j'^/2T^)e{T - Tc) where = 
l//3c. 



4.4- Density of zeros 

Now we calculate the density of zeros by using formula (Q. The density of zeros in each 
phase is given by 

f PI 

P/3(/3i,/32)= JV2n P2 . (48) 
[ SG 

In addition to this, zeros are distributed on phase boundaries where the derivative 
of the free energy shows a discontinuous change. On the P1-P2 boundary, that is, 
/3f + = /3c/2 and < |/3i|//3c < 1/2, the density of zeros is calculated as 



P/3(/3i, f^2) = 7TZ\7^ + 1^] {4>PiO{4>Pi - M + 4>P2e{(t)F2 - (49) 



= ^^m + (il-(ilm. (50) 

where 9 is the step function such that 6{x) = 1 when x is positive and when negative. 
In the same way, on the Pl-SG boundary \(3i \ + |/32| = /3c and 1/2 < |/3i|//3c < 1, 

P/3(/3i,/32) = ^mwi\-(3c+m). (51) 

Zn 

On the P2-SG boundary |/3i|//3c = 1/2 and 1/2 < |/32|//3c, the boundary distribution of 
zeros is shown to be zero. We conclude that the density of zeros in the REM is given by 

P/3(/3i,/32) = ^0iPj2 - |/3i|)^(/3? + P', - P!/2) 

+ ^ f^(/3i + /32 - P!me{pj2 - IAD 

+ + m - Pc)em - ^meip, - m. (52) 

Ztt 
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The result coincides with that in 
4-5. Discussions 

We have re-derived the density of zeros by using the rephca method, which allows us to 
clarify the following properties of the distribution: 

• Figures of the distribution of zeros can be viewed as phase diagrams in the complex 
plane. There are two types of distributions for the density of zeros: two-dimensional 
distribution within a phase and one-dimensional one at a phase boundary. 

• When (Jab = 0, the two configurations {ia} and {i'^} become independent of each 
other. Then, the average partition function decouples as (JT6ll and (p is written as a 
sum (f) = </)o(/3) + 0o(/3*)- In this case, the density of state is shown to be zero. 

• When Qab 7^ 0, there are correlations between two configurations {ia} and {i'a}- 
Then, two-dimensional distributions of zeros can be obtained. Since the correlation 
occurs only when the ensemble average is taken, one can say that the two- 
dimensional distribution is specific to random systems. In the case of the REM, 
two configurations are completely identical: {ia} = {i'a} in the P2 and SG phases. 

• However, qab 7^ is not the sufficient condition for two-dimensional distributions. 
In the REM, a two-dimensional distribution appears in the P2 phase and not in the 
SG one. Our result indicates that two-dimensional distributions are not specific to 
the SG phase. In the SG phase of the REM, the system freezes to its ground state 
and the entropy goes to zero. In such a frozen phase, zeros do not appear. 

• The two-dimensional distribution of zeros in the P2 phase is located far from the real 
physical axis and near the imaginary axis. On the imaginary axis, the Boltzmann 
factor is given by e~^^^^ and can take arbitrary values. Therefore, it is a natural 
result that zeros appear around the imaginary axis. It is interesting to see that 
the zeros can appear only when (32 > (3c/ on the imaginary axis. Since e~*^^^ 
is nothing but the time evolution operator, the result may be related to dynamical 
properties of the system. 

• The one-dimensional distribution of zeros on the Pl-SG boundary meets with the 
real axis at the angle of 7r/4, and the density of zeros is zero at the phase transition 
point on the real axis. These properties are characteristics of second-order phase 
transitions [27] . We discuss these properties closely in section [61 

• To understand the role of each zero, we calculate the specific heat from the obtained 
density of zeros. Using ([8]), we can write the specific heat as 



which allows us to write the function as a sum of contributions from each phase 
and boundary. The result is plotted in figure |H The singularity at the transition 
point comes from the zeros in the P1-P2 boundary. At large temperature, the main 
contribution comes from the zeros in the P2 phase. In the SG phase lower than 




(53) 



Replica analysis of partition-function zeros in SG models 



12 



the transition point, the contributions from the boundaries give a negative specific 
heat. All contributions cancel out to give the zero specific heat. As we decrease the 
temperature, zeros far from the real axis in the P2 phase become more important. 

• In section 14.21 we have defined the entropy for complex temperature by using 
the Legendre transformation. It is not obvious whether we should impose the 
nonnegative condition to the entropy since the function has not been defined as 
the number of states. In fact, the condition ( 1371) in the PI phase has not been 
used to determine the phase boundary. The results in the P2 and SG phases 
may be understood the number of states since the entropy is independent of the 
imaginary part of the temperature. Therefore, although the definition of the entropy 
in section 14.21 looks like a natural one, we need more discussion for the physical 
meaning of the entropy. 

5. GREM and the full RSB limit 

In contrast to a naive expectation that the two-dimensional distribution of zeros is 
specific to the SG phase, our result shows that zeros are not in the SG phase but in a 
P phase. The area-distributed domain is not on the real axis and cannot be directly 
related to the SG phase transition. 

The REM is considered the simplest SG model since the SG phase is described by 
the IRSB solution and the order parameter only takes either or 1. We want to discuss 
more complicated situations keeping the advantage of the REM that the analytical 
calculation is possible. For this purpose, we employ the GREM [231 1211 [251 126]. This 
model is a simple generalization of the REM and allows us to analyze higher step RSB 
solutions easily. We consider the distribution of zeros in the GREM and discuss possible 
distributions at the full RSB limit. 

5.1. Model 

The GREM is defined by a hierarchical structure of energy levels. Each energy level 
is expressed by the sum of K random numbers. To the i^th level of hierarchy with 
1 < I' < K, we assign random variables eu{l),e^{2), ■ ■ ■ ,ei,{M^). The number of 
variables M^, is given by 

= (54) 

where each is an integer with 1 < < 2^ satisfying 

Mk = {ai---aK)^ = 2^. (55) 

For the uth level of hierarchy, we generate random numbers with Gaussian distribution 
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where > satisfies 

K 

E«- = l- (57) 
From the random numbers generated as above, we construct 2^-random numbers as 

^^ = Ee.(L(^-l)M,/2^J+l), (58) 

where i = 1, 2, ■ ■ ■ , 2^ and [a; J is the floor function which indicates the largest integer 
not exceeding x. 

The GREM is deflned as a system with energy levels ( l58l) . By choosing parameters 
in an appropriate way, we can flnd multiple phase transitions. We deflne 

where z/ = 1,2, If Ti > T2 > ■■■ > Tk, phase transitions occur at these 

temperatures. 

In the replica analysis of this model, the SG order parameter q^^J is deflned at each 
hierarchy u [26] . At T > Ti, all the parameters are given by the P solution q^i^ = 5ab- 
At T = Tj,, q'^fj^ turns into the IRSB solution and the degrees of freedom in the z^th 
hierarchy fall into the ground state. For example, when we consider the K = 2 case, 
the system goes into the P-P, SG-P and SG-SG phases as the temperature is decreased. 
It is understood from the complexity analysis that the SG-P phase is identifled as the 
IRSB state and the SG-SG phase as the two-step RSB one. 



5. 2. Distribution of zeros 



It is a straightforward task to apply the calculation of zeros in the previous section to 
the GREM. We consider two hierarchy {K = 2) system and choose parameters so that 
the relation Ti > T2 is satisfled. We calculate = [In |2'|]/A^ in each phase. Following 
the calculation in |26l, we flnd 



ln2 + i(/32-/3|)j2 



I In «i + f /32 J2 + In «2 + fWf - Pl)J^ 
iln2 + i/32j2 

/3iJVailnai + lna2 + f{/3f - [5l)J^ 
/3iJVailnai + |lna2 + f(3lJ'^ 



PI - 
P2 - 
P2- 
SG 



PI 
PI 
P2 
-PI 



(60) 



SG-P2 
SG-SG 



/3i J [^Jai In tti + A/a2 In 0:2) 

which is easily understood as a straightforward extension of the result in the previous 
section to the two-hierarchy case. Then, comparing these solutions, we obtain the 
phase diagram and distribution of zeros in flgure [51 Zeros are distributed in the shaded 
regions and on the bold lines of the flgure. We have two phase transitions at T = Ti 2 
and correspondingly, two lines of zeros touch the real axis. Still, the shaded regions are 
far from the real axis and are not related to the SG transition. 
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Figure 5. Phase diagram 
and distribution of zeros of the 
GREM with two hierarchies in the 
complex- /3 plane. The plot is only 
in the first quadrant. Zeros are dis- 
tributed in the shaded regions and 
on the bold lines. 



Figure 6. A continuous limit 



is: cx) of the GREM. Critical 



temperatures change continuously 
between Ti = 1/(3^ and Tk = 0. 



An interesting result is obtained when we consider the continuum hmit of hierarchy. 
Then, the critical temperatures become continuous between Ti and and the SG phase 
can be understood as the full step RSB state. We can find a similar, but different, 
behavior as the Sherrington-Kirkpatrick model [2S] which also has the full RSB SG 
phase. In the GREM, the continuum limit is not unique and depends on the choice of 
parameters [21] . We show a typical example with Tx = in figure O We see that the line 
distributions at phase boundaries are accumulated to become area distributions. In the 
P2 phase of the figure, the original area distributions and accumulated lines overlap with 
each other and cannot be distinguished. The area distribution in the SG phase comes 
in contact with the real axis, which is consistent with the picture that phase transitions 
occur continuously in the full RSB SG phase. Actually, a similar numerical result has 
been obtained for a system on the Bethe lattice [12]. Our result in the GREM shows 
that distributions of zeros in SG systems with a full RSB state are more complicated 
than those in the REM. 

6. Ferromagnetic interaction and first order transition 

In the previous sections, we have studied the REM. This model is described by a 
random distribution of energy levels and the spin degrees of freedom do not appear 
at all. Therefore, in order to extend our method to other spin models, it is necessary to 
treat the spin Hamiltonian. To achieve this, we treat the p-spin model in (fT2|) . At the 
same time, to motivate the use of the spin representation, we treat the positive regular 
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interaction Jq whose effect is incorporated to the probabihty distribution of Ji-^^-.-ip as 

2' 



NP 



-1 



NP- 



exp 



^l^2■■■^p 



NP-^ 



(61) 



At large Jq, the system is in the ferromagnetic (F) phase which is characterized by the 
magnetization m = /N where ( ) denotes the thermal average. At the REM limit 

p — )■ oo, the model can be solved exactly, and the phase diagram in the Jq — T plane is 
obtained as in figure [3 The phase boundary between the SG and F phases is given by 
-^0 = -^oc = Jvhi2 and that between the P and F phases by 
P 1 

T = , (62) 

2 Jo-^Ji-jnn2 

The transition between the P and F phases is of first order and the entropy shows a 
discontinuous change. Therefore, by treating this model we can study both the spin 
representation and the first-order phase transition in our replica formulation. 
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Figure 7. Phase diagram of the 
REM with the F interaction Jq. 
The P-F transition is of first order 
and the others of second order. 



Figure 8. Phase diagram and 
distribution of zeros of the REM 
at Jo > '^Oc = J-\/ln2. The dashed 
lines are the phase boundaries at 
Jo < Joc- 



Following the standard method to calculate [Z"] [H Ej [3l [18] , we can treat in 
a similar way. Since the replica space is doubled to the complex conjugate space, we have 
two kinds of spin variables {S^} and {5"-"}. Correspondingly, the SG order parameters 
are defined as g.^ ~ EAStS'^)/N, q'^, ~ E^{S'/'S'^) /N and g,, ~ EAStS'^)/N. 
Since we treat the F phase, we also have the magnetization nia ~ J2i{Si)/N and 
m'^ ~ J2i{Si"')/N. By using these variables, we can write 

a>b 
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- N{p - Jo Y^i^aY - N{p - l)/3Vo EK)' 

a a 

+ AT In Tr exp Ip'j' ^labSaSb + E ^abS'aS'b 

y a>b a>b 

a,b a a J j 

where variables with the hat symbol are related to those without hat as 

qab = liqabY-'. (64) 

The order parameters are obtained from the saddle-point equations. At p — )■ oo, the 
equations are easily solved since the order parameters take either or 1. The PI, P2 and 
SG phases without the magnetization are calculated in the same way as the previous 
case and ( H5ll is obtained. In the F phase given by nia = rn'^ = 1, all spins are aligned 
to the same direction and we find a simple expression 

= (3iJo. (65) 

Comparing in each phase, we can write the phase diagram and the distribution 
of zeros in the complex /3 plane. When Jq < Jqc, the F phase solution is irrelevant and 
we have the phase diagram in figure [3l The result at Jq > -^oc is shown in figure M The 
phase boundary between the PI and F phases is given by 



f32J = v//3iV2-4/3iJo + 41n2 (66) 
and that between the P2 and F phases by 



(3iJ' = Jo-^J^-J'ln2. (67) 

The zeros are distributed in the P2 phase and on all the phase boundaries. 

Comparing figures |3] and [H we can see the difference between the second- and 
first-order phase transitions. The general consideration shows that there is a relation 
between the order of the phase transition and the distributions of zeros on the real 
axis [27]. In the second-order transition, the line distribution meets with the real axis 
with the angle of 7r/4, and the density of zeros goes to zero on the real axis. In contrast, 
for the first-order transition, the angle between the line and the real axis is 7r/2 and the 
density is finite on the real axis. Our results support these expectations. On the real 
axis, the density of zeros is given by 

J2 



P/3(/3i,/32 = 0) = - J2in2 9{Jo - Jv^) 

X 6 (^(3iJ^ - 2(Jo - pi - J2 In 2)^ . (68) 

It is also known from the general argument that the density at the transition point 
is directly related to the dis continuous e nergy change Ae. In the present case, it is 
obtained as Ae = ep — ep = J Jq — J'^ln.2. We see that the same factor appears in 
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7. Lee- Yang zeros of the REM 

Up to here we have discussed the Fisher zeros in the REM. In this section, we treat 
the Lee- Yang zeros in the same model. While the analysis of systems with complex 
temperature clarifies the system's thermodynamic properties such as the entropy and 
specific heat, that with complex magnetic field is useful to understand the magnetic 
properties such as the magnetization and susceptibility. 

The analysis of the Fisher zeros using the replica method goes along the same line 
as that of the Lee- Yang zeros in the previous sections. To apply the magnetic field, we 
consider two possible patterns: longitudinal and transverse magnetic fields. The density 
of zeros in the complex longitudinal field has been studied in [211 122]. The method used 
there is very similar to the original analysis of the Fisher zeros in the REM [19]. We 
re-derive the result by using the replica method. Also, we treat the system with a 
transverse field, which allows us to study the quantum fiuctuation effect. 



7.1. Longitudinal magnetic field 

We use the energy representation of the REM in ( ITO|) and (fTTj) . Then, the magnetic 
field h is incorporated to the partition function as 

N{M) 

Z = Y^e^hM j2 e-/3^K^^), (69) 

M 1=1 

where M represents the "magnetization" and takes 

Mk = -{N-2k), (70) 

with k = 0,1, ■ ■ ■ , N. For a given Mk, the number of the sum N{Mk) is equal to 
N\/k\{N — k)\. The partition function is written as 

z = e-^^'^Y.y' E e-^"""' (71) 

fc=0 i=l 

where y = e^^^ and E^i is a random variable taken from (fTTj) . Thus, the partition 
function is a polynomial of the A^th degree in y and the density of zeros is obtained 
from ([6]). 

The REM with real h has been treated in [Ml [17] and the phase diagram is 
depicted in figure [9l The P and SG phases are separated by the critical temperature 
T^{h) = 1/I3,{h) defined by 

^^Ml.e-.V,||m^. (72, 

As in the previous cases, this transition is also caused by the entropy crisis which means 
that the entropy goes to zero in the SG phase. 
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Figure 10. Phase diagram and 
Figure 9. Phase diagram of the distribution of zeros of the REM 

REM in the h-T plane. The P- in the complex-/i plane. The upper 

SG transition is of second order. gg^j^e is for Tc < T < 2Tc and 

lower for 2Tc < T. For T < T^, all 
domain belongs to the SG phase. 



When h is a complex number, the PI, P2 and SG phases appear in the phase 
diagram. is written as 

|^|2n _ g-27Vn/3/.i ^ ^ (^^ELo (y*)ELo KU^.}) 

{ka,ia} {k'aA} 

( N N{Mk) ^ 

xexp Yl M{ka,ia}) + n,,{{k'^,z'J))EkA, (73) 



k=0 i=l 



where 



n 



nkii{ka,ia}) = Y5kak5i,i, (74) 

a=l 

n 

nk{{ka}) = Y.^kak- (75) 

a=l 

Each ka takes 0, 1, ■ ■ ■ , and ia takes 1, 2, ■ ■ ■ , N{Mk^). Taking the average, we obtain 

[|^]2n ^ g-2iVn/3/xi ^ ^ (^)Er.o (/)EL, Ka/^^}) 

X exp I J2 ((lab + q'ab + '^Qab) > , (76) 

where the order parameters are defined as 

Qab = ^kakt^iaii, (77) 

Q'ab = hakMi',^ C^8) 

gab = Sk.kAaii- (79) 

As we mentioned above, we have three phases. 
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• PI. 

The simplest solution is given by qab = q'ab = ^ab and qah = 0. We have for 
= [ln|Z|]/A^ 

0P1 = -Ph, + ^ + \ ln(l + e^'^) + \ ln(l + e^^^*). (80) 

This is a simple extension of the usual P phase result for real h to the case of 
complex h. 

• P2. 

The second simplest case is specific to the complex parameters: Qab = q'ab = Qab = 
Sab- Then, we have the P2 phase result 

<f>P2 = -f3h + ^ + l ln(l + e'^'^). (81) 

This is valid at /3 < /3c(^i) = /3c(^i)/2 where the entropy is nonnegative. 

• SG. 

When the temperature is lower than l//3c(/ii), the system freezes to the SG phase. 
The order parameters are given by ( !39l) and f j4T|) . Then, 

0SG = -^h + ^-MMl. + in(i + e^^^e^^)'^^). (82) 

Comparing these results, we can draw the phase diagram in the complex-/;, plane 
as figure [T0][§|. Zeros can appear in the P2 and SG phases, and the P1-P2 and Pl-SG 
boundaries. The result is completely the same as that in the previous study obtained 
by a different method [21] . 

It is interesting to see that the density of zeros in the SG phase is nonzero as 

tanh^(2/3c(/ii)/i 
1 + 2(/i2/J2)(l -tanh2(2;ge(/ii)/ii)) 
in contrast to the case of the Fisher zeros where the zeros do not appear in the SG 
phase. In the present case, although the entropy is zero in the SG phase, the critical 
temperature depends on the magnetic field and the value of the free energy depends 
sensitively on the field. As a result, the density of zeros gives a finite value. But we note 
that the function (!83|) decays exponentially in hi and takes considerably small values, 
which mean that zeros very rarely exist in the SG phase. 

Another interesting point is that the P2 phase appears at T > 2Tc = l//3c(0). This 
property is the same as that of the Fisher zeros (see figures [3l [5] and [6]) . At present, the 
reason why 2Tc plays a special role is not clear. It may be interesting to see whether 
this property is specific to the REM or not. 

§ The phase diagram should in principle be drawn in the complex-?/ plane, though the translation 
between y and h is immediate. Then, the /i-plane is restricted to |2/3/i2| < tt. 
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7.2. Transverse magnetic field 

Then, we consider the REM in a transverse field [29] . The Hamiltonian is given by 

N 



H 



0", 



n<---<'T; 



i=l 



B4) 



where (T^ and are the Pauh matrices. Since two terms in the right-hand side of 
the equation do not commute with each other, the system must be treated quantum 
mechanically. At the REM limit p — oo, this model can be solved exactly and the phase 
diagram is given by figure [TTl We have two kinds of P phases: the classical P (CP) 
and quantum P (QP) phases. The CP phase is the same as the P phase in the classical 
REM and spins are directed along the z-direction randomly. On the other hand, in the 
QP phase, spins are aligned with the transverse-field direction. 
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Figure 11. Phase diagram of 
the REM in a transverse magnetic 
field r. The CP-QP and SG-QP 
transitions are of first order. 
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Figure 12. Phase diagram and 
distribution of zeros of the REM in 
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The average replicated partition function is written in terms of order parameters 
Xa(r, t') ~ Ei{SUr)SUr'))/N and q^r, r') ~ E^{SUr)S'Ar')) /N where a ^ b. Since 
we are treating the quantum system, the trace is defined by the imaginary time path 
integral with some measure [30], and the spin variable depends on the time r. The 
order parameter Xa with Qab < Xa < 1 characterizes quantum effects since Xa = 1 in the 
classical system. 

When r is real, the model is solved by the IRSB and static ansatz. We expect that 
the same ansatz can be applied to the case of complex F. Proceeding in a similar way 
as the previous cases, we have 

[\Zn = exp _ _ (1 _ m)q^) 



(p-l)(x'^-(l-m)g'^) 

{p-l)(x^-{l-m)qP)+N\nTTeA, (85) 



Replica analysis of partition-function zeros in SG models 21 

Nnl3*^J^ 
4 

Nn\f3\^J^ 
2 

^ a=\ B Wb / 

p2 72 n p2 j2 rijm' / m' ' ^ 

+ ^ (x' - 9') E(^:'^)' + E E 

^ a=l B \ae_B 

+ - 1) E s:s': + ^1 E E ^: E 

a=l ^ B \aeB J \a£B J 

n n ^ 

+/3r E + /3r* E ^1^^ , (se) 

a=l a=l J 

where S*" = f^ dr S'^{t)//3 and the order parameters with the hat symbol are defined in 
the same way as in ( !64|) . Possible saddle-point solutions are given as follows. 

• CP. 

X = x' = 1 and the other parameters are set to zero. In this case the quantum 
effect is irrelevant, and we have the CP result 

0cp = ^ + ln2. (87) 

From the nonnegative condition of the entropy, this solution is shown to be effective 
at T > Tc = J/2yhi2. 

• QP. 

When the transverse field is strong enough so that the effect of interaction is 
negligible, there is not any order in the 2;-direction. Then, all the order parameters 
are neglected, and we obtain 

0QP = ^ ln(e^^ + e~^^) + ^ ln(e^r* + e"'^^*). (88) 

• SG. 

At low temperatures the system freezes to the ground state. In this case, x = 0' = 1, 
and we have the standard IRSB result 



0SG = /3JVln2. (89) 

We see that any phase which is specific to the complex parameters does not 
appear. For example, we can examine the CP2 phase defined by x = x' = X = 1 
and q = q' = q = 0. Then, we have 

B^P 1 

0CP2 = ^ + 2ln2. (90) 

This solution is valid at T > 2Tc = J/ Vln 2. In this region, (/>CP2 is smaller than 0cp- 
Therefore, we can conclude that the CP2 phase is irrelevant. In the same way, we can 
also show that the "QP2" phase does not appear. 
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The phase diagram is shown in figure [T21 We do not have any two-dimensional 
distribution of zeros and zeros are distributed on the phase boundary. The phase 
boundaries between CP/SG and QP phases are determined by 

cosh(2/?ri) + cos(2/3r2) = ^c^'^cp/sg, (91) 

Their transitions are of first order, and we see that the general argument, the lines of 
zeros are vertical to the real axis, also holds in this quantum case. 

8. Conclusions 

We have discussed the partition-function zeros in random energy models by using a new 
replica method. We re-derived previously obtained results and obtained several new 
results on the distribution of zeros of random systems. Our method is very general and 
systematic and can in principle be applied to any spin model. 

We find that the two-dimensional distribution of zeros is characterized by the 
order parameter qab which connects configurations between the original and complex 
conjugated spaces. In the REM with complex temperature, the area-distributed phase 
is paramagnetic rather than a SG phase. It does not appear at real (3, but the physical 
quantities of the real system are affected by this phase any way as we have demonstrated 
for the specific heat. 

We also find differences between first- and second-order transitions. We showed that 
the general argument discussed in |27] holds even for random systems. Since our result 
is not a rigorous proof of the statement, it would be interesting to find a counterexample 
in other models. 

In contrast to previous numerical analyses of finite-dimensional SG systems, we find 
in the REM that the area distribution of zeros does not approach the real axis. This is 
because the REM is an oversimplified SG model. In the REM, the SG phase is described 
by the IRSB ansatz and the entropy goes to zero in the SG phase. These properties do 
not hold in general. In fact, the analysis of the GREM at the continuum limit shows 
a more complicated behavior. Therefore, the next thing to do will be to apply our 
method to other mean-field models such as the Sherrington-Kirkpatrick model [28] and 
the spherical model [3T| 132] . The formulation is the same as in the present case, and we 
can utilize the expression ( l63l) for general p-spin systems. The difference arises when we 
solve the saddle-point equation. The SG order parameters take values between and 
1, and possibly complex values. Such a situation does not exist for the REM, and we 
need a careful analysis, which will be done in a future work. 

Another point to be discussed is to see how the thermodynamic limit is obtained 
by taking the size of the system to infinity. We want to know how the zeros approach 
the real axis of the complex plane. Our analysis is based on the saddle-point method 
and is valid only in the thermodynamic limit. Therefore, we need to take a different 
way. Recently, the spherical model has been solved for arbitrary system size [33] • If 
the method developed there is also useful for systems with complex parameters, we can 
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study finite size effects for the distribution of zeros. It will be an interesting task to 
solve that problem in future works. 
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